Supersaturated dispersions of rod-like viruses with added attraction 
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The kinetics of isotropic-nematic (I-N) and nematic-isotropic (N-I) phase transitions in disper- 
sions of rod-like /rf-viruses are studied. Concentration quenches were applied using pressure jumps 
in combination with polarization microscopy, birefringence and turbidity measurements. The full 
biphasic region could be accessed, resulting in the construction of a first experimental analogue of 
the bifurcation diagram. The N-I spinodal points for dispersions of rods with varying concentrations 
of depletion agents (dextran) were obtained from orientation quenches, using cessation of shear flow 
in combination with small angle light scattering. We found that the location of the N-I spinodal 
point is independent of the attraction, which was confirmed by theoretical calculations. Surprisingly, 
the experiments showed that also the absolute induction time, the critical nucleus and the growth 
rate are insensitive of the attraction, when the concentration is scaled to the distance to the phase 
boundaries. 

PACS numbers: 82.70.Dd, 64.70.Md, 82.60.Lf, 83.80.Xz 



I. INTRODUCTION 



A long standing issue in the physics of fluids is the be- 
havior of the homogeneous fluid close to the point where 
it becomes unstable and phase separates, i.e. the spin- 
odal point ip s [l[ . Before the spinodal point the fluid will 
be meta-stablc or supersaturated, which means that the 
fluid will only undergo a phase transition when fluctua- 
tions in the concentration are sufficiently high to over- 
come a certain nucleation barrier. Thus the meta-stable 
region is characterized by the induction time T in< i for 
phase separation to set in. rind goes to infinity enter- 
ing the meta-stable region from the stable region at the 
binodal point (fib , i.e. 1/T m d — > at (pf,, while Ti n d — > 
approaching ip s . For molecular fluids it is very difficult 
determine the spinodal point, because the tiniest impu- 
rity will lower the nucleation barrier and the phase sep- 
aration is very fast. The binodal point is easier to ac- 
cess, since it is given by the final phase separated state. 
Colloidal systems have proven to be very suitable for this 
type of fundamental studies. The main reason is that the 
interactions between the colloids can be tailored, while 
the size of the colloids slows down the kinetics as com- 
pared to fluids permitting direct visualization One 
way of tailoring the interaction between colloids is to add 
non-adsorbing polymers to the system. Polymers induce 
attractive interaction between the colloids, due to the de- 
pletion of the polymer between the colloids [1] . The range 
and the strength of the attractive potential is controlled 
by the polymer size and concentration, respectively. Due 
to attractions colloid-polymer mixtures typically show a 
gas-liquid-like phase transition [J] . Despite of the advan- 
tages of colloids, it is also for this class of systems difficult 
to access the meta-stable and unstable region in a con- 
trolled way. Firstly, the spinodal and binodal line are 
located very close to each other, while the energy barrier 



for phase separation to take place is low. Secondly, sys- 
tems might be arrested in the metastable state@. The 
challenge is to bring the system in a meta-stable state 
while maintaining the system homogeneous. 

Where concentration is the only order parameter of 
interest for gas-liquid phase separating colloidal spheres, 
dispersions of colloidal rods have two order parameters 
that characterize the phase behavior which are strongly 
coupled: particle concentration and orientational order. 
These systems exhibit an isotropic-nematic (I-N) phase 
transition, where the system gains positional entropy at 
the cost of the orientational entropy The location of 
the transition can be derived from Onsager theory for 
slender hard rods Q. The phase behavior around the 
I-N transition is characterized by two branches 0: an 
isotropic branch with zero orientational order and a ne- 
matic branch with a finite orientational order, as depicted 
in Fig. [TJ The I-N binodal point ip{p will be first encoun- 
tered when following the isotropic branch by increasing 
the concentration. For ip > tp b the system will be meta- 
stablc, or supercooled, with respect to fluctuations in the 
orientation towards an aligned state. After some typical 
induction time T mt \ fluctuations will be sufficient to over- 
come the nucleation barrier and isolated nematic droplets 
(factoids) will grow in a isotropic background. For even 

higher concentrations, i.e. for ip > (p^p the system be- 
comes unstable to fluctuations, such that each fluctuation 
in the orientation of the rods will result in a continuous 
growth of the nematic phase out of the isotropic phase. 
Likewise, the system will have a N-I binodal ip b and 

spinodal ipi N ^ point when following the nematic branch 
by decreasing the concentration. In this case the system 
is superheated because fluctuations towards a lower or- 
dering drive the phase separation. Beyond the binodal 
point, i.e. for ip < (pi isotropic nuclei (atactoids) will 
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form after some induction time r; n d, while beyond the 
spinodal point, i.e. for ip < (pi , each fluctuation in the 
orientation of the rods initiates spinodal phase separation 
and an isotropic phase of disordered rods continuously 
grows out of the nematic phase. 

The difference in concentration between the I-N and 
N-I binodal points of a hard rod liquid is only about 
10 % 0. Hence the spinodal and binodal points are lo- 
cated very close to each other. The density difference 
between the binodal points increases when the rods are 
made attractive: the I-N binodal shifts to lower concen- 
trations while the N-I binodal shifts to higher concen- 
trations. This has been shown experimentally for vari- 
ous rod-polymer mixtures [§, 

[MEL 111- It is, however, 
not obvious that the location of the spinodal points are 
equally affected by adding polymer. As a consequence 
it is not known to what extend the metastable and un- 
stable regions are affected by the addition of polymer. 
Even for hard rods the location of the spinodal points, 
i.e. the open symbols in Fig. [TJ has never been experi- 
mentally confirmed. The goal of this paper is to locate 
the spinodal points with respect to the binodal points 
over a range of attractions by probing the kinetics of the 
phase separation process. This goal requires a well de- 
fined time t = at which the system is quenched from an 
initially stable state into a meta-stable or unstable state, 
in order to determine the induction time Ti n d. Taking 
advantage of the two order parameters that characterize 
rod dispersions also two types of quenches can be made: 
a quench in the orientation and in concentration. 

An orientation quench can be performed by first apply- 
ing an external field to a phase separated system some- 
where in the biphasic region, thus preparing a full ne- 
matic phase, as indicated by the dashed line in Fig. Q] 
The system is quenched when at t = the external field is 
switched off, which renders the system either unstable or 
meta-stable depending on the concentration. Tang and 
Fraden used the diamagnetic anisotropy of fd virus to in- 
duce FN phase transitions with a high magnetic field (HI 
and demonstrated nicely the existence of an unstable re- 
gion. Here they quenched, however, always to some fi- 
nite field strength. Similarly, we used in an earlier paper 
shear flow to prepare a stable nematic phase [l4| . Ces- 
sation of flow at t = renders the nematic phase meta- 
stable or unstable depending on the concentration, see 
Fig. [T] Small Angle Light Scattering (SALS) was then 
used to probe the formed biphasic structures and deter- 
mine the spinodal as the concentration where structure 
formation sets in immediately after the quench. In this 
paper we again rely on this technique on mixtures of fd 
and dextran, but as compared to the earlier experiments 
the sensitivity is improved such that measurements could 
be performed also when density differences between the 
phases were small as is the case at low polymer concen- 
trations. 

The disadvantage of the orientation quench is that only 
the nematic-isotropic transition is probed. In order to 
access also the I-N transition one needs to make a con- 



centration quench. Such a quench can be made by rig- 
orous stirring a phase separated system and probe it im- 
mediately after stirring, as was done for dispersions of 
boehmite rods [l5j . 

Apart from the practical problems this imposes on the 
experiment, the results could also be biased by resid- 
ual alignment in the sample after mixing [l6| . Quenches 
were also initiated by polymerizing short actin fila- 
ments [13 ■ Both methods could evidence the distinction 
between nucleation-and-growth and spinodal decompo- 
sition by the morphology of the phase separated struc- 
tures, without pinpointing the actual location of the spin- 
odal points. Nuclcation- growth mechanisms and spin- 
odal structures have also been observed in computer sim- 
ulations [TH, [l^, H3, HH • In this paper we perform pres- 
sure quenches from 1 bar up to 1000 bar and vise versa. 
Given the compressibility of water, this corresponds with 
instantaneous concentration quenches of up to 5 % [22| . 
as indicated by the solid arrows in Fig. [1] We probe 
changes using polarization microscopy, birefringence and 
turbidity measurements. We could determine both the 
isotropic to nematic (I-N) and nematic to isotropic (N- 
I) spinodal, since full nematic phase could be induced, 
starting with a full isotropic phase. Thus we construct 
a first experimental analogue of the bifurcation diagram 
plotted in Fig. [T] Since with pressure quenches only a 
small concentration range can be accessed, we relied on 
cessation of shear flow to study the attractive rods with 
added dextran. 

To supplement the experiments we have used Scaled 
Particle Theory (SPT) approach to predict the phase be- 
havior of colloidal rods for a range of polymer concen- 
trations, including the I-N and I-N spinodal lines. The 
experimental data are qualitatively compared with this 
theory. In Sec. |TT] we introduce the SPT and present 
its results in the form of two phase diagrams. In Sec. 
IIIII wc introduce our experimental techniques and sample 
preparation. The effect of the pressure and orientation 
quenches are given in Sec. UlTBl and ImCl respectively, 
resulting in the determination of the spinodal and bin- 
odal points in Sec. IIVCI and growth rates in Sec. IIVDI 



II. THEORY 

The phase diagram of a rod-polymer mixture can be 
predicted from free- volume theory, as elaborated in detail 
in Refs. [H, [24|. The free energy per particle of a system 
of N hard spherocylindrical rods with length L and di- 
ameter D in a volume V in osmotic equilibrium with a 
reservoir of ideal, non-adsorbing polymer with a volume 
fraction ifp takes the following form: 

(i) 

in terms of the thermal energy (3 1 = ksT, rod as- 
pect ratio 7 = L/D 1 and polymer-colloid size ratio 
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Figure 1: Bifurcation diagram of the nematic order parameter 
S for hard rods (ip^ = 0) with L/D = 133 corresponding to 
the free energy Eq. ([1]). Open circles indicate spinodal points, 
while the filled circles indicate binodal points. The arrows 
indicate concentration quenches that are made to render the 
system in a supersaturated state. The dashed line indicates 
possible locations of the system when inducing a full nematic 
phase applying an external field. 



q = 2R g /D (with R g the polymer radius of gyration). 
The density variable y = ip/(l — ip) is related to the rod 
packing fraction ip = (n/4)LD 2 N/V. The reference part 
ifip = corresponds to a system of hard rods and stems 
from Scaled Particle Theory (SPT). The orientational en- 
tropic contribution in Eq. |T]) is defined as 



*[/] 



du/(u)ln[47r/(u)] 



(2) 



where the unspecified distribution /(u), describes the 
probability of rods with orientational unit vector u, nor- 
malized over all possible orientations via J du/(u) = 1. 
The coefficients P and Q pertain to the shape (i.e. aspect 
ratio) of the rods: 



P[f] = 3 

Qif] = 



3(7 ^If 
37-I 

127(27- 1) 



r[f] 

, 127(7 -I) 2 
(3 7 -l) 2 + (3 7 -I) 2 



[/] (3) 



The quantity r[/] represents the following double orien- 
tational average: 



r lf] = ~[[ dudu'f(u)f(u')\u x u' 



(4) 



The last term in Eq. (JTJ) accounts for the depletion con- 
tribution. It depends on the free volume fraction a, ex- 
pressing the average fraction of the system volume avail- 
able to the polymer at a given rod packing fraction ip. 
An explicit expression follows from SPT: 

a([f], <p) = (!-*>) cxp (-Ay - B[f]y 2 - C[f]y 3 ) (5) 



with coefficients A,B,C given explicitly in [2J]. Since 
the reservoir polymer concentration tpR is proportional 
to the depth of the minimum of the attractive depletion 
potential, it serves as a measure for the strength of attrac- 
tion between the rods. As the polymers are treated as an 
ideal gas, the polymer volume fraction ip po i y in the sys- 
tem simply follows from multiplying the reservoir value 
(fip with the fraction of available free volume «([/], ip). 

The SPT coefficients P, Q and free volume fraction 
depend implicitly on the unspecified orientational distri- 
bution /(u) . In the isotropic state, all orientations are 
equally probable so that / = l/4n, a = and r = 1. 
In the nematic state, it will be a non-uniform distribu- 
tion peaked along some nematic director n. An accurate 
variational form for / has been proposed by Onsager Q 
which takes the following form: 



/(*) 



/icosh(recos( 
47r sinh k 



(6) 



with < 6 < 7r the polar angle between u and the 
nematic director n (cos 9 = u • n) and k > a varia- 
tional order parameter (note that k = leads back to 
the isotropic constant / = l/47r). With the use of an ex- 
plicit trial function, the orientational averages associated 
with Eq. ((4| and Eq. ([2|) become analytically tractable 
[fiLHH- The results for the orientational averages are: 



<t(k) = ln(KCothK) — 1 



arctan(sinh k) 
sinh k 



and 



2J 2 (2«) 

t[k) = ; 5 — — 1 

2 sinh k 



>0 (7) 



(8) 



with l2{x) a modified Bcsscl function. 

The equilibrium value for At is found by a minimization 
of the total free energy which leads to the stationarity 
condition: 



dF 
dn 







(9) 



for any given rod packing fraction and attraction strength 
ifip. The nematic order parameter associated with the 
equilibrium value for k is found from: 



S = 



duV 2 (u-n)f(u) 



= 1 - 



3 coth k 



(10) 



where S = in the isotropic and < S < 1 in the 
nematic phase. The solution of Eq. ([9]) for hard rods 
(ifp = 0) is given in Fig. [1] showing two branches where 
the stationary solutions correspond to a local minimum 
of the free energy 0]. The spinodal points marks the 
transition between a stable and an unstable solution of 
Eq. ©. 

However, this (second order) transition is preempted 
by a first order phase transition. Thus a discontinuity 
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both in concentration and in the ordering of the system 
will occur. The co-existence of two phases requires that 
the osmotic pressure II and chemical potential /i of the 
isotropic phase, with volume fraction (p\!\ an d the ne- 
matic phase, with volume fraction (rf , are equal: 

i%«) = i%W) (ii) 
= M N) ) (12) 

The binodal points can thus be found using the ther- 
modynamic relations /i = 
combination with Eq. [TJ 

The phase diagram for a rod-polymer mixture is given 
in Fig. [2^ and shows the characteristic widening of the 
biphasic gap as the amount of polymer is increased. The 
location of the spinodal points, however, appears much 
less affected by the depletion attraction. This is reflected 
more clearly if we plot the spinodal curves in terms of the 
fraction < (p nem < 1 of ncmatic phase formed, rather 
than the overall rod packing fraction ip. Applying the 
lever rule, we may relate <p ne m to ip via: 



§7jv,t and n 



BF 
dV N,T 



Pu 



(N) 
<Pb 



4 I] 



(13) 



with p^ N ^ the binodal rod packing fractions corre- 
sponding to the coexisting isotropic and nematic phases. 
Fig. [2Jd shows that the NI spinodal instability occurs if 
the overall rod concentration corresponds to the nematic 
phase occupying about 20 % of the system volume. This 
result is virtually independent of the strength of the de- 
pletion attraction as long as the amount of added poly- 
mer is not too large. 



III. EXPERIMENTAL 
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Figure 2: (a) Phase diagram for a rod-polymer mixture with 
L/D = 133 and colloid-polymer size ratio q = 2R g /D — 5.4. 
Plotted in terms of the rod packing fraction tp TO d and polymer 
volume fraction ip po iy in the system. Coexisting isotropic and 
nematic phases are connected by tie lines, with the nematic 
phase having a higher ^ r od- (b) Location of the isotropic- 
nematic spinodals in terms of the fraction of nematic phase 
jp ncm [Eq. (|13[) ] plotted versus the system polymer concentra- 
tion yp iy on the vertical axis. 



A. Sample 



Fri-virus suspensions were used in a 20 mM Tris buffer 
with 100 mM NaCl at a pH of 8.2. The virus is a long 
and thin rod like particle (length 880 nm long, width 
6.6 nm, persistence length 2.2fim). Attractions between 
the rods are varied through depletion by addition of dex- 
tran (480 kd, Pharmacosmos). A small amount of FITC- 
labeled dextran was added to be able to determine the 
dextran concentration spectroscopically. See Ref. [26| 
for the labeling procedure of dextran. The samples were 
prepared as follows: First, a homogeneous fd- virus sus- 
pension of 21.1 mg/ml fd virus with dextran is allowed 
to macroscopically phase separate into an isotropic and 
nematic phase. The two phases were then separated into 
two different vials and the dextran and fd-virus concen- 
trations were determined spectroscopically. Three differ- 
ent dextran concentration were used in this study given 



an initial virus concentration of 21.1 mg/ml: 6 mg/ml 
(low), 13 mg/ml (middle) and 20 mg/ml (high). The re- 
sulting phase diagram is shown in Fig. [3] This phase 
behavior differs somewhat from previous published re- 
sults for the same system 0, [l4|. The deviation might 
be due to different polydispersity of the dextran which 
can drastically change the interaction 27[ and thus the 
phase behaviour [28j. By combining different volumes of 

the isotropic, tp^ , and the nematic, tp^ , phase from the 
initially phase separated sample we can prepare any con- 
centration along one tie-line, with a concentration tp ncm 
relative to the phase boundaries as expressed in Eq. 1131 
The concentration of dextran and fd-virus for each sam- 
ple was checked after every new mixing. For pressure 
experiments a sample very close to the isotropic-nematic 
spinodal has been prepared at the same ionic strength, 
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Figure 3: Phase diagram of the I-N transition of [dextran] 
vs. [fd] at an ionic strength of 110 mM. The solid symbols 
are the binodal points as determined spectroscopically after 
phase separation. The open symbols are the spinodal points 
as determined after shear rate quenches, see below. The thin 
lines connecting the binodal points are the tie lines. 



3 hours. Optical birefringence measurements were per- 
formed on a home made setup consisting of a argon ion 
laser line at 514. 5nm, a beam expander, two polarizers 
and a Coherent Fieldmaster-GS power and energy me- 
ter. The detector was operated by a self written Lab- 
VIEW application. The pressure quenches were applied 
with the same protocol as described above, starting al- 
ways at 1 bar. The detector registered the intensity of 
the through-going beam with a time resolution of 0.4 s 
over a total time of 1 hour. To acquire more information 
about kinetics of the nematic-isotropic transition, mea- 
surements of the forward transmission were performed. 
When a phase transition takes place, the sample becomes 
turbid which is indicated by a decrease in the transmis- 
sion. Turbidity was probed with a Hc-Nc laser and a Co- 
herent Fieldmastcr-GS power and energy meter. Pres- 
sure releases were performed as described above. The 
transmission was detected and registered as in the bire- 
fringence measurements with the help of the Lab VIEW 
software. 



but without dextran. 



B. Microscopy, Birefringence, and turbidity at 
high pressure 

For all pressure experiments we used a small container 
sealed by a vitron ring, which contains the sample while 
it allows for pressurizing via holes in the brass support 
ring. For microscopy and birefringence measurements 
the container was placed in the polarization microscopy 
cell, where the polarization was maintained. The cell 
was mounted in a specially designed cell holders for mi- 
croscopy or birefringence. For turbidity measurements a 
SANS cell was used which has longer optical pathway, to 
increase the sensitivity of the experiments, but which has 
windows that scramble the polarization. The detailed de- 
scription of the cell and the preparation procedure can 
be found in Ref. [2<|. 

Polarization microscopy was performed on a Zeiss Ax- 
ioplan microscope with a 10x/0.30 Plan - NEOFLUAR 
objective. The cell was placed between crossed polariz- 
ers in order to detect birefringence corresponding to the 
change in the order parameter. Starting at 1 bar we 
applied different pressures up to 1000 bar with steps of 
200 bar steps. Pressurizing was performed with a rate 
of 100 bar/s. Pressure releases were performed from a 
sample in the nematic phase that had been equilibrated 
at a pressure of 1000 bar for about 1 h. Consequently, 
we decreased the pressure starting from this equilibrated 
sample inducing the nematic-isotropic phase transitions. 
After changing the pressure, sequences of images were 
taken using a CarlZeiss Axiocam CCD colour camera and 
Zeiss acquisition software. The sequences of images were 
taken with different time resolution: in the beginning ev- 
ery 5 s and in the end every minute over a total time of 



C. Small Angle Light Scattering (SALS) under 
shear 

A homcbuild optical couette shear cell combined with 
a SALS-setup was used [30| . The shear cell consisted of 
a rotating inner cylinder with a diameter of 43 mm and a 
static outer cylinder with a diameter of 47 mm resulting 
in a gap width of 2 mm. The inner and outer cylin- 
ders were both made of optical grade glass. A 15 mW 
HeNc Laser (Melles Griot) operating at a wave length of 
632 nm was used. To ensure that the laser beam went 
through only one gap, it was directed through the cen- 
ter of the rotational axis of the inner cylinder. In the 
rotating cylinder, the beam was directed along the ra- 
dial direction with a prism. Scattered intensities were 
projected with a lens directly on to the chip of a Peltier 
cooled 12-bit CCD camera, with 582x782 pixels (Prince- 
ton Instruments, microMAX). The scattering angles on 
the chip were calibrated by placing a known grid (PAT 
13 Heptagon) in the scattering volume. 

The fd-virus solution was always pre-sheared at 
100 s _1 and quenched to zero shear rate at t = 0, at 
which time the registration of the SALS patterns started 
with a rate of about two frames per second. Immedi- 
ately after cessation of flow the rods will on average be 
oriented along the flow direction. Thus initially the sys- 
tem will be in a homogenous flow-induced nematic state 
with a well defined 'director' n. The effect of the pre 
shearing was checked for a number of samples by reduc- 
ing the pre-shearing to 50 s . No difference on the result 
could be found and we therefore kept the protocol of a 
pre-shearing rate of 100 s" 1 . 
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1 bar to: 400 bar 800 bar 1000 bar 



Time 
s 

5 min 
15 min 

40 min 
3 h 




Figure 4: Sequence of polarization microscopy image after 
an increase in pressure, starting at 1 bar. The final pres- 
sure is a measure of the concentration. Bright regions with 
a higher orientational order parameter appear at about 15 
minutes after a quench to 400 bar indicative of a nucleation- 
growth mechanism (middle left). For a final pressure of 1000 
bar phase separation sets in immediately (top right), indica- 
tive of spinodal decomposition, while after 40 minutes a full 
nematic phase is formed (bottom right). 



IV. RESULTS 

A. Concentration quenches using pressure 

Sequences of micrographs taken after pressure 
quenches for different time delays are gathered in Fig. 
2] and The starting pressure for Fig. 2] was 1 bar 
so that initially the system was in the isotropic phase. 
The varied depth of the quenches allowed for exploration 
of different regions on the phase diagram. For the 200 
bar pressure quench the images stay dark throughout 
the whole experiment (data not shown). For a quench 
to400 bar the total intensity only starts to increase af- 
ter about 5 minutes. At longer times brighter regions, 
indicative of a finite order parameter, are visible due to 
nucleation-and-growth. The possible presence of factoids 
could not be observed, because a lOx objective was used. 
An induction time Ti n d for phase separation to set in is 
characteristic for the meta-stable states. For a quench 



1000 bar to: 600 bar 400 bar 200 bar 




Figure 5: As Fig. |4]but now starting at fOOO bar. Nucleation- 
and-growth events can be seen for a final pressure of 600 bar 
after 15 minutes (middle left). For a final pressure of 200 bar 
phase separation sets in immediately (top right), indicative of 
spinodal decomposition, but also the total intensity is much 
reduced, indicative of a lower ordering. After 1 hour a full 
isotropic phase is formed (bottom right). 



to 800 bar an almost instant change in the intensity is 
observed along with homogeneous structure formation. 
This observation suggests that the I-N spinodal is located 
in the immediate vicinity of this applied pressure. In 
the last sequence corresponding to the deepest quench of 
1000 bar one can see the full transition from isotropic to 
nematic phase through spinodal decomposition. First, 
the increase in intensity started immediately after the 
quench, i.e. Ti n d = 0. Second, the initial early stages 
of the transition exhibited morphology characteristic for 
spinodal decomposition - interconnected, labyrinth-like 
structures spanning through the whole sample. Figure 
[5] shows a thresholded picture taken at 50 seconds after 
this quench. Typical spinodal decomposition morphol- 
ogy is evident. The last two micrographs of the 1000 
bar quench show that the system develops into a fully 
nematic phase indicated by large homogenous regions of 
the same colour that only reorient with time. Thus the 
FN spinodal tp"p is located in the very proximity of the 
N-I binodal <p[ N) . 

Fig. [5] depicts the nematic-isotropic phase transition 
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Figure 6: Filtered and binarized micrographs of fd virus after 
the pressure jump from to 1000 bar taken at 50 seconds 
after the quench. Growing interconnected structures are a 
clear proof of spinodal decomposition taking place. 



after the pressure quenches from 1000 bar to the indi- 
cated pressures. The difference in the morphology for 
the shallowest quench from 1000 to 800 bar are subtle 
and changes are only visible 3 h after quenching (data 
not shown). For the quench from 1000 to 600 bar the 
changes are more pronounced. After about 20 minutes 
grainy structures are formed, characteristic for a phase 
separation via nucleation and growth. Such structures 
can still be seen in the last image of the sequence though 
the overall intensity decreased substantially. The image 
taken two minutes after the quench from 1000 to 400 bar 
shows that the overall intensity has decreased, indicative 
of a lower ordering, while already some biphasic dark 
and white structures seem to be visible. The first image 
after the deepest quench, from 1000 to 200 bar, shows 
that phase separation sets in immediately (top right), 
indicative of spinodal decomposition, while the total in- 
tensity is much reduced. After 1 hour a full isotropic 
phase is formed (bottom right) showing that a full phase 
transition took place. A more conclusive location of the 
spinodal and binodal points can be obtained with turbid- 
ity and birefringence measurements as presented in Sec. 

Since with 1000 bar we approach the limit of the max- 
imum applicable pressure, we cannot use this technique 
to study dispersions of rod-polymers mixtures. For such 
dispersions the width of the biphasic region (see Fig. [3]) 
cannot anymore be bridged. For this reason we use shear 
flow to induce a flow stabilized nematic phase followed 
by cessation of shear flow. 



B. Orientation quenches using cessation of shear 
flow 



The development of the scattering pattern after a 
quench in the shear rate can be seen in Fig. [7J for three 
different if ncnl : 0.20, 0.62 and 0.82. The time at which 
the first detectable scattering structure appears increases 
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Figure 7: Development of the scattering pattern for 
<^ncm =0.20, 0.62 and 0.82, after a shear rate quench from 100 
to s _1 for the high dextran concentration. For if nC m = 0.20 
spinodal decomposition immediately sets in as can be con- 
cluded from the faint ring at 31 s (upper right), while for 
V^ncm = 0.82 nucleation of structure can only be observed 
after about 150 s (middle left). The indicution time for struc- 
ture formation is indicated by the dashed line. 



with increasing ^ ncm , as indicated by the dashed line 
in Fig. [7] For the highest concentration this is the 
case only after about 150 seconds. For this concentra- 
tion it is difficult to identify specific features. For the 
lower two concentrations the scatter pattern occurs as 
two slightly bend lines perpendicular to the director (x). 
In the early stage these lines scatter weakly but with 
time the intensity increases, the peak sharpens and the 
position moves to smaller wave vectors. The scattering 
patterns are isotropic in the initial stage but become in- 
creasingly more asymmetric with time. This indicates 
that the formed objects are increasingly more elongated 
and oriented along the director. To quantify the develop- 




time [s] 



Figure 8: Development of the scattering along the flow direc- 
tion for </?ncm = 0.20 (a) and ^ ncm = 0.62 (b) as deduced from 
the scattering patterns in Fig. [JJ The total time is three and 
a half minutes in (a) and six minutes in (b). 



merit of the formed biphasic structure we extract inten- 
sity profiles along and perpendicular to the director. The 
time development of the scattering along the director is 
shown in Fig. [5) Here the formation of the structure peak 
and both the movement to lower q-values and the inten- 
sity increase of this peak are clearly seen for the two low- 
est concentrations, <p n em, where it should be mentioned 
that for the lowest ip nem the peak develops immediately 
after cessation of flow, indicative of spinodal decompo- 
sition. At the highest concentration the structure peak 
can only be distinguished at longer times, characteristic 
for nucleation-and-growth. Fig. [5]and similar plots taken 
for other concentrations will be used in the following sec- 
tions as a base to deduce the induction time for structure 
formation. 



C. Spinodal and binodal points 

Figures Hl5l and [71 all show that after the quench phase 
separation sets in immediately or after some induction 
time Ti n d, depending on the depth of the quench. As ex- 
plained in the introduction, in order to locate the spin- 
odal and binodal points we have to find at what con- 
centration ri n d for the formation of nematic (for the I-N 
transition) or isotropic structures (for the N-I transition) 
goes either to zero or to infinity. 

To determine T; n d for quenches in the orientation we 
plot the intensity of the peak in the scattering pattern, 
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Figure 9: (a) The intensity at the peak of the scattering pat- 
tern, J((?max), after cessation of shear flow for j; ncm = 0.30 
(□), tpnem = 0.62 (o) and tp n em = 0.82 (a) at the high dex- 
tran concentration to study the N-I transition, (b) Turbidity 
(solid symbols) and total intensity of the polarization micro- 
graphs (open symbols) after a pressure drop (cubes: 400 bar, 
bullets : 200 bar) to study the N-I transition. The total in- 
tensity of the polarization micrographs after the initial decay, 
/total j i s a l so indicated, (c) Birefringence intensity after an 
increase in pressure to study the I-N transition. The lines in 
the plots indicate the extrapolation to determine the induc- 
tion time after a shear rate or pressure quench. 



see Fig. [5J as a function of time. Fig. [5£i shows this time 
development for three different j3 nem . The induction Ti n d 
is now obtained by extrapolating the linear intensity in- 
crease to zero intensity. This procedure was repeated 
for different attraction strength, i.e. different dextran 
concentrations. The result is shown in Fig. [TUb . where 
respectively r; n d and l/rind arc plotted vs. </? nein for the 
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three attractions. The N-I spinodal tfs and binodal 
P om ts ar e determined by extrapolating T- m ^ and 
1/ T ind to zero, respectively. Interestingly the curves of 
Ti n d and 1 /ri n d f° r the different attractions overlap if we 
scale the concentration relative to the phase boundaries. 
As a consequence the location of ipi is independent of 
the attraction and is found to be ip nem — 0.25. More 
striking even is the observation that also the absolute in- 
duction times are only effected by the relative distance 
from the phase boundary tp ncm and not the absolute fd 
concentration or the attraction, i.e. the dextran con- 
centration. Of course the phase boundaries can also be 
determined by measuring the concentration of the two 
separated phases or visual observation, considering that 
no scattering pattern should occur for the homogeneous 
phases. The deviation between the different methods is 
less than 5%. This means that the high concentration 
phase boundaries can be determined in three different 
ways and the low concentration phase boundary in two 
ways. 

The spinodal decomposition is characterized not only 
by the fact that it immediately sets in, but also by the 
morphology of the formed structure, which should be bi- 
continuous. It has been shown that the transition in the 
morphology is quite smooth [Til , [l5| . This is also exem- 
plified by Fig. [7] and Fig. [5b for ip ncnl = 0.62. A clear 
ring structure is observed, exemplary for spinodal struc- 
tures, but also for this concentration an induction time 
is observed. Thus the phase separation has already spin- 
odal characteristics, but clearly the system is still in the 
meta-stable region. Therefore the only reliable way to 
determine the spinodal point is by using the induction 
time. 

The kinetics of the phase separation for samples with- 
out polymer was studied using the pressure quench. The 
reason is that with this technique we access both the N- 
I and the I-N transitions. As described in section Ull Bl 
the location of these points could be estimated from mi- 
croscopy experiments. The interpretation of the micro- 
graphs is, however, not straightforward. Quenches start- 
ing from an initial nematic phase at 1000 bar show that 
the overall intensity decreases as well as that structures 
are formed. To separate the two effects we additionally 
performed turbidity measurements, since turbidity is a 
measure of the biphasic structure that is created in the 
sample during phase separation. In Fig. the responses 
of the turbidity and the total integrated intensity Itotai 
of the polarization micrographs are plotted for quenches 
to 400 and 200 bar. For the quench to 400 bar clearly 
an induction time is observed in the turbidity, while the 
birefringence, which dominates the intensity of the po- 
larization microscopy, shows a fast and a slow decay. 
Since the flow induced nematic phase has a higher or- 
der parameter than the metastable nematic branch (see 
Fig. [J) the nematic phase will first relax to this branch, 
which explains the initial fast decay of the birefringence 
to ^to™J ch . Turbidity is not sensitive for this process, 
because structure formation is not involved in this pro- 
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Figure 10: Induction times (filled symbols) and inverse in- 
duction times (open symbols) starting with an initial nematic 
phase (a,b) to determine and c J c>i JV ' ) and with an initial 

isotropic phase (c) to determine ip^ and ipi 1 ^ ; (a) as a func- 
tion of ifinem for three dextran concentrations (high (□), mid- 
dle (o) and low (a)) after pre-shearing; (b,c) as a function of 
the final pressure after an initial pressure of 1000 bar and 1 
bar, respectively. The lines are a guide to the eye. 



cess. Hence we can conclude that phase separation sets 
in only after an induction time of about 650 seconds. For 
the quench to 200 bar changes in birefringence and tur- 
bidity set in immediately so Tj n d = 0. Since we know 
from microscopy that the final stage for 200 bar is fully 
isotropic, the down turn in the turbidity after half an 
hour can be interpreted as the disappearance of biphasic 
structure. Before reaching the fully isotropic phase (low 
turbidity) the system undergoes a phase transition (high 
turbidity). When pressurizing an initially isotropic sam- 
ple then the increase of birefringence can only be caused 
by the formation of the nematic phase. Thus the induc- 
tion times for the formation of the nematic phase can be 
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Figure 11: The experimental equivalent of the bifurcation di- 
agram. The pressures corresponding to <p[ N ' > and ip~P are de- 
termined from Fig. [9b and c as the pressures where rind — ► 
(solid symbols). The pressures corresponding to and ip^ 
are given by the pressures where 1/rjnd — > (open symbols). 
The total intensity of the polarization micrographs after the 
initial decay, itoSu > ( see Fig El 3 ) is taken as a measure of the 
orientational order parameter (stars). The line is a guide to 
the eye. 



obtained from the birefringence responses as plotted in 
Fig. He. 

The induction times and the reverse induction times 
are plotted in Fig. [TTjb and c for the turbidity and bire- 
fringence, respectively. The pressures corresponding to 



and tpi 1 ^ are determined from these two figures as 



the pressures where t\ 



Ind — > 0, while the pressures cor- 
responding to tpfr and ip^ are given by the pressures 
where l/rjnd — » 0. We can now construct the experi- 
mental equivalent, Fig. [Til of the theoretical bifurcation 
diagram, Fig. [JJ using pressure as a measure of concen- 
tration and the total intensity of the micrographs after 
the initial decay as a measure of the orientational order 
parameter (Fig. [5Jd). Note that the error bar in the 
pressure is determined from the uncertainty in the ex- 
trapolation to T; n d — > and l/rj n d — ► 0. Another source 
of error could also be the exact concentration of the rods, 
since the data were taken on different batches. This er- 
ror would show up as a shift of the entire curve that is 
obtained from the used batch. This could explain why 
ipi is somewhat smaller than , which is in principle 
not possible. 



D. Growth rates 

The average size of the structures that are formed dur- 
ing the phase separation can be deduced from the loca- 
tion of the peak of the scattered intensity. The resulting 
size parallel and perpendicular to the director arc plotted 
in Fig. [T^] as a function of time. Surprisingly the size of 
the formed structures as well as the rate with which they 
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Figure 12: Time development of the structure size parallel 
(solid symbols) and perpendicular (open symbols) to the di- 
rector, for three dextran concentrations at a scaled volume 
fraction of ip n em = 0.30. Symbols as in Fig. 110b . 
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Figure 13: The growth rate vs. (pnem parallel (open symbols) 
and perpendicular (solid symbols) to the director for three 
dextran concentrations. Symbols as in Fig. HOfa , 



grow do not depend on the attraction. The structure 
formed in the early stage is isotropic, while the struc- 
ture growth is anisotropic: the structures are growing 
faster along the director so that the structure becomes 
anisotropic with time but the anisotropy in the struc- 
ture does not seem to increase beyond an aspect ratio of 
2. This can also be appreciated from Fig. [T31 where the 
extracted growth rates are plotted against (p ncul for all at- 
tractions in. The only factor that affects the growth rate 
is the distance from the phase boundaries. As expected 
the growth rates go to zero at </4 while the growth 
rate seems to go to a plateau value around i^ nem = 0.25, 
which is, as we have seen in the last paragraph, the spin- 
odal point ipi^ . 

To access the very early changes, below 25 s, we look at 
the logarithm of the intensity ln(7) with time, at different 
the wave vectors q. In this way we probe the length scale 
of the density fluctuations which has initially the highest 
probability to grow. This investigation is only done along 
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Figure 14: ln(Int)/dt vs. g for a) three different tp ngm : 
Vncm = 0.20 (□), (f ncm = 0.30 (o) and tp ncm — 0.70 (A) 
for the middle dextran concentration and b) three dextran 
concentrations at j3 nHn = 0.2. Symbols as in Fig. HOfa , 



the director since for the neutral direction the signal to 
noise ratio is too low. In Fig. Wfin d\n(I)/dt is plotted 
as a function of q for three different ip ne m for sample B. 
d\n(I)/dt peaks at the same q vector for all <p n cm, i-e. 
independent if it is in the unstable region or far in the 
meta-stable region. The length scale over which the sys- 
tem grows fastest turns out to be of the same size as the 
objects which have just formed at the induction time, 
i.e. the critical nuclei. For the nucleation-growth regime 
this means that in the first time window new nuclei arc 
formed continuously with a well defined size, which start 
growing at the induction time. In Fig. [Tib d\n(I)/dt vs. 
q is plotted for the three different attractions at a volume 
fraction of tp nem = 0.2, i.e. for a spinodal decomposing 
sample. Again the maximum q of the growth is the same 
for all attractions, as was found for the induction time. 
On the other hand the amplitudes arc different. This am- 
plitude different is due to the change in contrast differ- 
ent between the isotropic and nematic. With increasing 
dextran concentration the width of the biphasic region 
increases and thus the density and contrast between the 
two phases. If we now combine this observation with the 
observation that both the growth rate of the objects and 
the induction time do not change with added attraction 
it is evident that what by eye seems to be a faster phase 
separation with more dextran added only is a optical ef- 



fect due to the increasing in contrast. 



V. DISCUSSION 

The main goal of our investigations was to study the 
extent of the supercooled or superheated regime for dis- 
persions of attractive rods and how it depends on the 
strength of the attraction. This dependence is given in 
Fig. [3j The location of the N-I spinodal point almost 
does not change in the range of attractions studied here. 
This holds both for the absolute concentration of fd virus 
(Fig. [3J as well as the fraction of the coexisting nematic 
phase </9 nem (Fig- 1 10b ). at least within the experimen- 
tal error. These observations are confirmed by theory, 
where it is found that both spinodal lines as a function 
of the absolute rod concentration are insensitive of the 
attraction (Fig. [2^). The N-I spinodal line has a weak 
dependence on the attraction when plotted as a function 

of (Fig. [2d). 

Using ip ncm as the scaled concentration we observed 
that not only the location of the spinodal is insensitive 
for the attraction, but also the size of the critical nuclei 
that are formed in the meta stable region, Fig. [T21 the 
induction after which they start to grow, Fig. 110b . and 
the rate with which they grow, Fig. 1121 This is surprising 
since in classical nuclcation theory these parameters de- 
pend on the difference in the chemical potential between 
the homogeneous and demixed state and the interfacial 
tension between the two formed phases [l[ . These thermo- 
dynamic parameters are expected to be quite different, 
considering that the width of the biphasic region for the 
highest used dextran concentration has increased with 
an order of magnitude, see Fig. [31 Given the fact, how- 
ever, that the relative location between (p^ and tp s 
set the length and time scales, there must be another pa- 
rameter that determines the nucleation barrier. To gain 
theoretical understanding of this problem would require 
a full dynamical density functional approach, including 
spatial inhomogeneities 3l|, to access the evolution of 
the microstructure over time. Using Monte Carlo simula- 



tions Schilling et al. [321 ] did not find an increased orienta- 
tion correlation for polymer volume fraction comparable 
to those used in the experiments. If the initial stage 
of phase separation is dominated by collective rotational 
diffusion, then both experiments and the simulation (32| 
hint that depletion interactions do not lead to a stronger 
preference of rods to align. Simulations by Cuetos et 
al.[l9j confirm that the aspect ratio of critical nucleus 
is less than two, independent of supersaturation. The 
size of the critical nuclei is however in the order of one 
rod length, whereas our measurements indicate that the 
critical nucleus is about 7 rods long independent of the 
depth of the quench or the attraction strength. This is 
most clearly shown in Fig. 1 141 because this figure plots 
at what density wavelength the intensity growths fastest. 
Combination of simulations and theory does show that 
the anisotropy in the surface tension due to the planar 
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anchoring of rods at the interface plays an important role 
in the formation of a critical nucleus [2(3]. If this effect 
dominates the kinetics, it is understandable that attrac- 
tion between the rods is less important than for example 
the aspect ratio of the rod. Comparing experiments and 
simulations one should realize that not only the aspect 
ratio of the viruses is an order of magnitude bigger than 
those used in simulations, but also that we consider in 
Fig. [14] the N-I and not the I-N transition. 

The only parameter that is influenced by the attraction 
is the rate at which the gradient in the density grows, 
as plotted in Fig. 114b . which is due to the increasing 
width of the biphasic region. All samples in this plot are 
quenched into the unstable region and therefore undergo 
spinodal decomposition. The curve for the highest at- 
traction, i.e. the highest polymer concentration, shows 
a clear peak. This is typical for spinodal decom pos ition 
that is dominated by translational diffusion [33l. |34| . At 
the middle concentration there is still a clear peak, but 
it is not obvious wether the demixing rate goes to zero 
for q — > or not. For the lowest concentration of dex- 
tran the data are too noisy to draw any conclusions. Our 
results are somewhat different from earlier experiments 
performed at 13 mg/ml of dextran, which hinted that 
the growth rate does not go to zero for q — > [Til ]. 
Also the binodal lines found in this paper have sharper 
features as compared to those published in Ref. [§]. Dif- 
ferences in the poly-dispersity of the dextran that is used 
could explain these discrepancies. In the latter paper 
also a significantly large part of the phase diagram was 
covered. The reason for the somewhat limited range of 
attraction studied here is that higher concentrations of 
dextran would lead to a too high turbidity of the sample 
and multiple scattering, corrupting the reliability of the 
measurements. This effect could explain the fact that 
for the high polymer concentrations in Fig. 114b a fi- 
nite growth rate is found at high q values. Concerning 
the time dependence of the structure growth we found 
within experimental error a power law of around one, 
see the linear dependence in Fig. |12b. whereas theory 
predicts a lower power dependence [31[. Possibly we are 
restricted to the very initial stage of the phase separation 
process. The sizes of the coalescing structure formed at 
later times are too big, so that the scattered light hits 
the beam stop. This problem does not occur when using 
microscopy as in Ref. [13, Hf| ■ 

With the pressure quench we accessed both the N-I 
and the I-N transition. As for the N-I transition we ob- 
served that also tpi 1 ^ is located in the proximity of ■ 
With this we confirm for the first time the theoretical 
prediction done in Ref. 0]. What these experiments 
also show is that for the deepest quenches the initially 
homogeneous single phase (I or N) undergoes a local 
phase separation before it completely turns in to the new 
single phase. This mechanism confirms similar observa- 
tions in computer simulations [20| on initially isotropic 
hard sphcrocylinder-polymer mixtures. The fact that we 
can reach the full ncmatic state from the isotropic state 



with the pressure quench is striking, since it is known 
that the concentration difference between the two bin- 
odal points is 10 %, whereas with a pressure quench to 
1000 bar the water is compressed only 5 %. Here it is 
important to note that we are comparing the isotropic 
phase at 1 bar with the nematic phase at 1000 bar. As- 
suming that the width of the biphasic region does not 
change with pressure, it would mean that at high pres- 
sure the I-N transition sets in at lower concentrations. It 
is known that the location of the I-N transition is tem- 
perature dependent. This dependency could be linked 
to the temperature dependence of the flexibility of the fd 
virus [36| . Similarly, the features of fd virus could changes 
at high pressures. This change cannot be an irreversible 
process like denaturation, since the phase transitions are 
completely reversible. Further experiments are needed to 
show if is a specific feature like the flexibility of rod-like 
viruses that changes with pressure. Alternatively, more 
general features like the Debye double layer of charged 
colloids could be pressure dependent. 



VI. CONCLUSIONS 

We have studied the behavior of supersaturated dis- 
persions of rod-like viruses. Superheated nematic disper- 
sions were prepared by first applying a strong shear flow. 
Cessation of the shear flow at t = renders the nematic 
phase meta-stable or unstable depending on the concen- 
tration, see Fig. [1] We probed the structure formation 
using Small Angle Light Scattering. With the analysis 
of the scattering patterns we could access the induction 
time for structure formation r; n d, the size of the criti- 
cal nucleus and the growth rate. These parameters were 
measured over a broad range of attractions as induced 
by the addition of dextran as a depletion agent. We 
found that the N-I spinodal point, i.e. the concentration 
where r; n d — > 0, is independent of the attraction, which 
was confirmed by theoretical calculations. Interestingly, 
also the absolute induction time, the critical nucleus and 
the growth rate are insensitive of the attraction, when 
the concentration is scaled to the distance to the phase 
boundaries as given by v? n cm- This observation hints 
that concepts of classical nucleation theory are insuffi- 
cient to understand nucleation processes in anisotropic 
fluids. We also applied pressure quenches on dispersions 
of rods without added polymer, thus supercooling or su- 
perheating the system. The pressure quenches were deep 
enough to induce a complete phase transition from the 
isotropic to the ncmatic phase and vice versa, which takes 
place initially via phase separation. As a consequence 
both the N-I and I-N spinodal could be accessed. By 
a combination of polarization microscopy, birefringence 
and turbidity measurements we were able to construct a 
first experimental analogue of the bifurcation diagram of 
Kayser and Raveche Q. 
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